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Abstract 

Relying on the Hidden Local Symmetry (HLS) model equipped with a mechanism 
breaking the U(3)/SU(3)/SU(2) symmetries and generating a dynamical vector meson mix- 
ing, it has been shown that a global fit successfully describes the cross sections for the 
e + e~ — > 7r + 7r~, e + e~ — > (vr°/r/)7 and e + e~ — > 7r°7r + 7r~ annihilation channels. One 
extends this global fit in order to include also the dipion spectra from the r decay, taking 
into account all reported information on their statistical and systematic errors. A model ac- 
counting for lineshape distortions of the spectrum relative to p° is also examined when 
analyzing the r data behavior within the global fit framework. One shows that a successful 
account for e + e~ annihilation data and r spectra can be simultaneously reached. Then, 
issues related with non-perturbative hadronic contributions to the muon g — 2 are exam- 
ined in details. It is shown that all e + e~ data considered together allow for improved and 
motivated estimates for the a M (7r + 7r~), the 7r + 7r~ loop contribution to the muon g — 2 ; for 
instance, integrated between 0.630 and 0.958 GeV, we find a M (-7r + 7r - ) = 359.62 ± 1.62 
(in units of 10 -10 ), a 40% improvement of the current uncertainty. The effects of the 
various r samples in the context of a global fit procedure leads to conclude that different 
lineshape distortions are revealed by the ALEPH, BELLE and CLEO data samples. Re- 
lying on global fits to the data quoted above, one also provides motivated estimates of the 
7r + 7r~, 7r°7, 777 and 7r°7r + 7r~ contributions to a M up to 1 GeV with the smallest possible 
uncertainties. These estimates are based on various global fit configurations, each yielding 
a good probability. 
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1 Introduction 



It has been proved in [1J that the scope of the HLS model ESI, suitably broken flU, can be 
extended in order to include annihilation processes like e + e~ — > (tt° /r)^, e + e~ — > 
and decay spectra like 77/77' — ^ 7r + 7r~7, beside the e + e~ — > n + ix~ cross section. Actually, it 
also includes e + e~ — > K + K~ and e + e~ — > K°K° annihilations which have to be examined 
separately as they raise known specific problems [5|. Actually, most VMD physics up to the 
<p meson mass is covered by our extended model [Q]|, except for channels involving scalar 
mesons or channels where higher mass vector mesons could have a significant influence (61 
as, seemingly, e + e~ — > tun. However, all the e + e~ annihilations channels examined in HI are 
reasonably well described up to the <p mass region by the model presented in HI. 

The issue is now to examine how the Extended (HLS) Model performs while including 
other processes like the r ± — > 7r ± 7r°z/ decay which is in the scope of the HLS model. This 
leads us to report on the results of global fits using the existing r dipion spectra beside the e + e~ 
data extensively discussed in [1]. The same issue was partly addressecQin our former 01. The 
energy range of our model is limited approximately by the 4> meson mass ; meanwhile, as far 
as issues like the muon g — 2 value are concerned, this is an energy region where a reliable 
model can address some questions in a novel way. 

Indeed, an interesting outcome of such a global fit is the estimate it provides for various 
hadronic contributions to the muon g — 2. This covers the ix + n~ loop contribution, but also 
those from the (n fr/)^ and final states. The improvements following from using 

simultaneously r data and e + e~ data as well as the consequences of having a reliable global 
model describing all VMD physics up to the <p meson mass are interesting issues. Indeed, the 
underlying unified (HLS) framework of our model correlates several decay channels because 
of their common underlying physics and phenomenological studies indicate that these physics 
correlations are well accepted by the data [1]. Stated otherwise, one can examine in great details 
several consequences of accounting for r decays and e + e~ annihilations within a consistent 
framework. 

This turns also out to readdress the long-standing problem of the discrepancy between the 
BNL measurement Q of g — 2 and the predictions based on e + e~ annihilations and r spectra 
as reported in the literature j8l|9l[T0l[IIl[l2l[l3]|. A quite recent study |[l4l tends to lessen the 
disagreement between these two kinds of predictions, but not to resorb it. 

Our present study is based on all the e + e data sets used in [OQ and on the published r 
spectra. These are the dipion mass spectra collected some time ago by ALEPH [fT5l and CLEO 
[fT6l ; a valuable data set collected by the BELLE Collaboration, with a statistics of several 
million events, has been made recently available ifTTl . 

The paper is organized as follows. In Section [2l we describe the model for the dipion 
spectrum in the r 1 * 1 — > ^ ±r K Q u decay. The vector meson mixing produced by Isospin Breaking 
(IB) BUIllI is complemented in a simple manner with another IB mechanism allowing lineshape 
distortions of the p ± meson mass spectrum compared with p°. In this Section, we also list the 
r data sets and outline how they intervene in the global fit procedure. 

In Section[3l we first emphasize the correlation between pure IB shape distortion parameters 

'We considered together with e + e~ — > ir + ix~ cross sections, the pion form factor from ALEPH and the dipion 
spectrum lineshape from CLEO. We were, thus, less sensitive to issues related with the absolute scale of the r 
spectra. 
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and the absolute scale of r spectra. Then, the consistency of the dipion spectra from r decay 
- not affected by vector meson mixing - with all e + e~ annihilation channels is investigated. 
The behavior of each r data set - the ALEPH lfT5l . CLEO lfT6l samples and the recently issued 
BELLE IfTTl data sample - are examined under various kinds of fit conditions. It is shown that 
the ALEPH spectrum can be well described with simple and intuitive IB lineshape distortions 
compared to e + e _ , whereas this does not work well with BELLE and CLEO spectra. The best 
way to account for these is rather a rescaling of the absolute scale of their spectrum. We argue 
that this could point towards a more complicated IB lineshape distortion model than ours. One 
shows, nevertheless, that a satisfactory simultaneous account of all e + e" annihilation data and 
the available dipion spectra in the r decay can be reached, under quite reasonable conditions. 

In Section HI we focus on the non-perturbative hadronic contributions to the muon g — 2, 
especially the ix + n~ one. The results provided by the various e + e" — > n + ix~ data sets are 
examined and the effects of global fits involving the e + e~ — > (n° /rj) f y and e + e _ — > 
cross sections is shown. The effects of including the r dipion spectra within the fitted data sets 
is examined in full details. It is also shown that the KLOE data set |fT8l for e + e~ — > n + ix~ does 
not lead to some inconsistency but, rather, allows for improved uncertainties at the expense of 
a worse fit probability. 

We also conclude on the most likely value of several contributions to the muon g — 2 
following from a global fit to a large sample of e + e~ and r data. The uncertainties yielded look 
much improved with respect to usual. 

Finally, Section[5]is devoted to a summary of our conclusions. In particular, one emphasizes 
using the various r spectra in order to provide - or improve - theoretical predictions for the 
muon g — 2, taking into account the difficulty to model lineshape distortions in a way accepted 
simultaneously by the ALEPH, BELLE and CLEO data sets. 

In the present paper, which is the second part of a study started in [1], one does not discuss 
the properties or the results of the fits to the e + e~ data in isolation. These have been discussed 
at length in [Q]| ; we also refer the reader to that paper for the details of our model. All notations 
have been carefully chosen in order to ensure full consistency with flTJ . Finally, the present 
work supersedes and improves large parts of our former study flU. We also correct here for a 
(minor) computer code error which affected the treatment of the sample-to-sample correlated 
uncertainties in the data sets from |fT9l [201 I2T1 ; this is, indeed, important in order to provide 
reliable uncertainties to our g — 2 estimates. 

2 Including r ± — > tt ± 7t°u Data 

The difference between e + e" and r based estimates of the hadronic contribution to the 
muon g — 2 is an important issue. Indeed, accounting for isospin symmetry breaking effects, the 
t dipion spectral provide predictions for the hadronic contribution which makes the expected 
value of g — 2 close to its experimental measurement 0. Instead, all theoretical estimates based 
on e + e~ data deviate by more than 3 a. Comprehensive discussions of this issue can be found 
in jHHHQil an d more recently in lfT3l . Summaries can also be found in lfTTl[T2l . for instance. 
A quite recent reanalysis of this discrepancy [14J concludes to a smaller disagreement between 

2 Each normalized to the world average branching ratio Br(r — > ttttv), highly influenced by the ALEPH 
measurement |fT51l . 
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r and e + e~ based approaches (about 2 a) ; consequently, the newly proposed r based estimate 
moves farther from the BNL measurement. However, even if reduced, the mismatch between 
e + e _ and r based estimates of the hadronic contribution to g — 2 survives. 

It was shown in [4] that an appropriate account of isospin symmetry breaking (IB), in- 
cluding its effects on the (p, to, 0) mixing, certainly solves a part of the reported discrepancy 
between e + e~ and r spectra. However, the IB vector mixing defined there and recalled in [OQ 
does not exhaust all effects of IB. In this paper, we examine more deeply than in 01 the effects 
of IB shape distortions and their connection with absolute scale issues. In order to examine this 
kind of IB, one needs a data sample where the p° (e + e~ annihilation) and the p ± ( r decay) 
spectra are simultaneously present. 

The problem of the hadronic contributions to the muon g — 2 was not addressed in H|. 
This issue is examined here in a wider context by revisiting the consistency pattern of the 
t 1 * 1 — >• n ± ir i / data on one hand, and the much larger data set on the e + e~ — > e + e~ — > 

(n°/r])^ and e + e~ — > annihilation channels on the other hand. This is allowed by 

having extended the model presented in 0J in such a way that anomalous and non-anomalous 
channels are implemented within the unified framework presented in (TJ. 

Most part of the e + e~ — > ix + n~ data sets has been commented on in fl4) ; here, we only 
remind how the sample-to-sample correlated part of the systematic uncertainties should be 
treated, as this plays an important role in estimating the uncertainty on the muon g — 2. All 
other e + e _ annihilation channels have been considered in details in our recent [1]. Because 
of the poor probability of the best fit to the KLOE data |fT8l already commented on in 0]|, 
the corresponding data sample is not included systematically in the set of e + e _ — > n + n~ data 
samples considered ; however, its effects will be commented upon at the appropriate places. 
Finally, in order to fit the parameters of our p, u, mixing scheme [4], one still uses a subset 
of 9 radiative decay width data which have been taken from the latest issue of the Review of 
Particle Properties ||22| and are given explicitly in 0]|. 



2.1 The Model For r ± -> ^tt^u Decay 



with 



Our model for the pion form factor in r decay coincides exactly with the formulae given in 



FZ(s) 



2 ; 



P v 



DJs) 



fI-Il w ( S ) , fl = agfl 



(1) 



(2) 



m 



n 



where Hw{ s ) accounts for the loop corrections to the p^ — W T transition amplitude f T p and 
IIpp(s) is the p ± self-mass. Both loops are such that they vanish at s = 0. g denotes, as usual 
B2l[Il, the universal vector coupling and m 2 = ag 2 f% is the p meson mass squared as it occurs 
in the HLS Lagrangian. a is the standard HLS parameter expected close to 2 [El |4l [TJ| . 

Beside the mixing of vector mesons produced by breaking Isospin Symmetry, Reference flU 
examined the possibility of having a mass difference beween the neutral and charged p mesons. 
Here, we also allow for a mass squared difference between neutral and charged p mesons - 
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denoted resp. m? and m? + 5m? . Additionally, we also allow for a coupling difference of 
these mesons, resp. g and g' = g + 5g. The p ± — transition amplitude should be modified 
correspondingly [4], as will be reminded shortly. These two parameters correspond within our 
model to allowing mass and width differences between the charged and neutral p mesons, as 
commonly done in other studies lfT4ll23ll . 



2.1.1 The Pion Form Factor In the r Decay 



With the IB modifications just defined, the pion form factor has to be slightly modified 
compared with Eq. It can be written : 



FZ(s) 



(1 



2' 



^D p {s) 



(3) 



where g' = ag'/2 = a[g + 5g]/2. The other ingredients are modified, compared with Eqs. 
CQ) and ©, and become : 



D P {s) 



r;-nw(s) 



s — m 



5m? 



f; = f T +sr , */: 



5m? 



(4) 



9' 



9' 



where f T p = agf 2 is the p — W transition amplitude, D p (s) is the inverse p ± propagator and 
II' (s) is the charged p self-mass. With the 5/J term in the last Eq. ©, F£(0) = 1 is 
identically fulfilled. In dUl, we assumed 5g = 0. 

The (modified) F^(s) is the W — p transition amplitude with its loop corrections. In terms 
of the pion ^(s) and kaon £k(s) amputated loops, one has the following expressions : 



U w (s) 



7 pTVTT 



[Sp7T7r 



+ Pw(s) 



24 



£k{s) 



(5) 



P P {s) 



where z& = [Ik/ fir} 2 is the standard SU(3) breaking parameter in the BKY breaking scheme 
[[24ll25l . while Pw(s) and P p (s) are subtraction polynomials with real coefficients to be fixed 
by external conditions. 

One could look for a motivated way, like the BKY mechanism (|24|. able to generate this 
kind of IB distortion effects. The proposed modifications look, however, reasonable and cor- 
respond to the usual way of introducing mass and width differences in other studies. This 
mechanism will be referred to as IB shape distortion and, if numerically relevant, may comple- 
ment the IB vector mixing flUEQ. 

We have checked that one can safely identify £ n (s) and £r(s) - both being charged-neutral 
meson loops - occuring in these expressions with the amputated -k + tt~ and K + K~ loops ap- 
pearing in e + e~ annihilations flU. 
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In order to reduce the number of free parameters in the global fit procedure, we still identify 
(as in [4J) the subtraction polynomial for Hw(s) with those for its partner in e + e~ annihilation 
(see Section 5 in 0]|). On the other hand, as one can neglect pseudoscalar meson mass differ- 
ences in loop calculations, one also identifies the charged p self-mass IT' (s) with its neutral p 
partner - up to the 5g effect - as reminded in Section 5 of [0Q. 

Finally, in order to fit r data, one has to correct for specific isospin symmetry breaking ef- 
fects. For this purpose, short range |[26l (S EW = 1.0235) and long range flZZHMIH (G EM (s)) 
radiative corrections have to be considered. While comparing with experimental data, the quan- 
tity in Eq. © has to be modified to : 

f;( s ) [s ew g em (s)} 1/2 f;( s ) (6) 

As the standard HLS model does not go beyond the lowest lying vector meson nonet, we 
cannot fit the whole dipion r decay spectrum. We chose to stop around the <p mass HI as, 
up to this energy, higher mass vector mesons seem to have a very limited influence within the 
channels examined in [1J and in the present work. 



2.1.2 Useful Expressions For The Dipion Partial Width in the r Decay 

The dipion partial width in the decay of the r lepton can be written flU : 



647r 3 m 3 



ds 



\F n (s)\ 2 G (s) + 0(e 2 ) 



(7) 



with 



G a (s) 



4 K-s) 2 K + 2s) 
3 s 3 / 2 ^ 

J\s - (m n o + m 7r +) 2 ][s - (m n o - 



(8) 



where e = (m 2 — m 2 r+ )/m 2 r+ ~ —0.06455. The terms of order e 2 - which manifestly break 
Isospin Symmetry - are negligible. On the other hand, one obviously has : 



1 dY T , 



ds 



1 dN(s) 
N ds 



(9) 



where is the (integrated) nn partial width in the r decay ; l/NdN(s)/ds is the normal- 
ized spectrum of yields over the accessible dipion invariant mass ranged. While referring to r 
normalized spectra in the following, we always understand this quantity. 

Using Eqs. CD) and © together with the customary expression [|22l| for the the r — > ev T v e 
partial width, one can derive : 



2ml 



1 B vn 1 dN(s) 



\V ud \ 2 (m 2 T -s) 2 (m 2 +2s)/3 3 B e N ds 



(10) 



3 Of course, the total number of pion pairs is defined by N = J [dN(s) / ds] ds. 
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which is the standard expression used by experimentalists to reconstruct the pion form factor 
from experimental data [16, 17J. In this expression /3_ is the pion velocity in the dipion rest 
frame, B^ and B e are the branching ratios of the r decays to resp. tvtvu t and to ev T v e . 
Eq. © can also be written : 

1 dT nn (s) _ 1 dN(s) 

T T ds ™N ds ' 

where r r denotes the full r width. The theoretical expression for dT^/ds on the left-hand side 
is given by Eq. © and by (^(s)) 2 as following from Subsection 12.1.11 above ; the additional 
factors shown by Eq. © are be understood. Finally, the numerical value for T T - not accessible 
to our model - is derived from the measured lifetime [|22l| and B n7T is numerically provided by 
each experiment with various uncertainties. Eq. (TTTT) is the main tool in the present analysis. 

2.1.3 Absolute Normalization of the Dipion r Spectrum 

As clear from Eqs. © and ©, the absolute normalization of the theoretical dipion par- 
tial width spectrum is determined by the product G 2 F \V u d\ 2 S ev/G em{s) . Correspondingly, the 
absolute normalization of the experimental spectrum on the right-hand side of Eq.dTTb is deter- 
mined by the branching ratio B^. 

Less obvious analytically, but numerically important, is the role played by the universal 
vector coupling g and the transition amplitude /' T in providing the theoretical normalization 
of the dipion spectrum. Indeed, as a is found numerically close to 2, Eqs. §3§ and (HJ) show 
that the absolute magnitude of the dipion spectrum is proportional to the product squared of the 
p — W and pun amplitudes. Therefore, actually, non-vanishing 5g and 5m 2 influence both the 
lineshape and the absolute normalization of the r spectrum. 

Moreover, one cannot exclude some other mechanism breaking CVC by modifying essen- 
tially the absolute normalization of the r spectra ; therefore, a correction factor (1 + r] C vc) ma Y 
enter Eq. (TTTT) and can be fitted. Related with this, one should note a recent BaBar measurement 
about the t — jj, — e universality. BaBar report^] 113011 g^/g e = 1.0036 ± 0.0020 as expected, 
while g T j g^ = 0.985 ± 0.005 exhibits a 3 o departure from 1. If confirmed, this may indicate a 
possible CVC violation in the r sectoiE] affecting only the absolute scale of the r spectrum. On 
the other hand, an experimental bias on the B nn branching ratio cannot be excluded and could 
play in the same direction. 

Even if very close to standard approaches, the IB lineshape distortions have been introduced 
here in a very simplified manner. One can easily think of a more complicated structure of these 
than the one we inferred. 

2.2 Dealing With The Fitted Data Sets 

Besides the data sets provided^ by the ALEPH lfT5l and CLEO lfT6l Collaborations, a new 
sample has been recently made available by the BELLE Collaboration[17J. These are the r data 

4 We thanks W. M. Morse to have drawn our attention on this paper. 

5 This BaBar result contradicts a former measurement from ALEPH [15] which was consistent with lepton 
universality. 

6 As in our 0, we discard the OPAL data set ||3T1 . 
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sets which will be examined in conjunction with the whole set of e + e~ data samples already 
considered in |QQ|. We remind that a subset of 9 vector meson decay partial widths is also 
used, corresponding to decay modes not related with the annihilation data considered in our fit 
procedure ; these are numerically extracted from G2l . 

2.2.1 The r Input To The Fit Procedure 

In the present study, we submit to fit the experimental spectra as shown in the right-hand 
side of Eq. (TTTI) . In order to remain consistent, we use for each experiment its own published 
branching ratio measurement and not the world average branching ratio. 

As for the CLEO data, our input is their published spectrum lfT6ll for l/NdN(s)/ds nor- 
malized to their latest updated branching ratio measurement lT32l |22| . B W7r = (25.36 ± 0.44)%. 
This Collaboration also claims an uncertainty on the absolute energy scale lfi~6l l33l of about 
0.9 MeV. However, in our former analysis JH, no such uncertainty showed up significatively. 
Anticipating somewhat on our present analysis, we confirm its effective consistency with zero 
and, therefore, discard this freedom from now on. 

Concerning the ALEPH data, we use directly the last update of the B^/NdN(s)/ds spec- 
trum Q3]|. The corresponding branching fraction, B n7T = (25.471 ± 0.097 ± 0.085)%, is the 
most precise among the published measurements. The uncertainties will be added in quadrature 
(0.127%). 

For the BELLE data IfTTl , we have been provided [341 with all information concerning 
the pion form factor spectrurry, its covariance matrix for statistical errors and its systematics. 
The systematics have been added in quadrature to the statistical error covariance matrix. The 
BELLE l/NdN(s)/ds spectrum data are published as such [17] ; its error covariance matrix 
can be derived from the corresponding information provided for the pion form factor, using 
simple algebra. As stated above, we have submitted to fit the BELLE B 7T7r /NdN(s)/ds spec- 
trum normalized to the BELLE branching ratio 03 B n7T = (25.34 ± 0.39)%. 

The uncertainty provided by the branching ratio error is clearly a scale uncertainty and a 
bin-to-bin correlated error ; this should be treated as reminded in Section 6 of HI. This turns 
out to define the (partial) \ 2 for each of the ALEPH, BELLE and CLEO data sets by 0]| : 

xlxp = K 1 + ^Ex P )mi - /(si)][(l + Xexp^j - /(sj)]^" 1 + 

having defined, for each experiment, the measurements rrii as the central value for the branch- 
ing ratio times l/NdN(si)/ds and V being the full error covariance matrix. f(si) should be 
understood as the left-hand side of Eq. (fTTj) computed at the appropriate energy point. 

For each experiment, A is a scale parameter to be fitted and r] is the ratio of the branching 
ratio uncertainty to its central value. The second term in this x 2 is the standard way to account 
for a scale uncertainty. We have ijcleo = 1-74%, rj ALEPH = 0.51% and rj BELLE = 1.53%. 

With this input to the fit procedure, the ALEPH, BELLE and CLEO data sets are clearly 
treated on the same footing. As emphasized in H), if for some experiment the ratio Xfit/vExp 
is small enough (typically not greater than ~ 1 -4- 2), one can neglect this scale correction and 

7 Normalized to the world average branching ratio Br(r — > tt + it~ v), four times more precise than the BELLE 
own measurement, and slightly shifted. 



^Exp 
VExp 



(12) 
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use the standard x 2 expression in the minimization procedure, with the replacement =>■ 
Vij + rf E miTTij. Otherwise, one may consider we are faced with some missing variance and 
keep Eq. (fT2l) as it stands. 

In a previous study flU, we limited ourselves to considering only the r data points up to 
0.9 GeV in order to avoid at most effects of higher mass vector mesons. In the present work, 
however, preliminary studies using the BELLE and CLEO datc(§ samples lead us to push this 
upper energy limit up to 1.0 GeV. Indeed, as in our e + e fit studies JT), the influence of higher 
mass vector mesons seems negligible all along the energy region from threshold to 1.0 GeV 
- and even slightly above ; therefore, there is no physical ground to abstain from such an 
extension of the fitting range. 

2.2.2 Testing The r Spectrum Lineshapes 

The remarks presented in Sub section [2X3] explain why it is certainly appropriate to test the 
various r spectrum lineshapes independently from their absolute magnitudes. This can be done 
in two different ways. 

A first method turns out to normalize the data points to the sum of the data points 
covered by our fitted energy range (from threshold up to 1 GeV/c). Then, correspondingly, the 
model function function f(s) on the left-hand side of Eq. (fTTI) should be normalized to its 
integral over the fitted range. 

Another method, is simply to minimize the x 2 as defined by Eq. (PT21) . but amputated this 
from the (A/?]) 2 term which constrains the scale in accordance with the claimed experimental 
uncertainty. Indeed, in this way, the scale factor is allowed to vary freely within the global fit 
procedure. We checked that these two methods give similar results. 

2.2.3 Dealing With The Uncertainties In e + e Data Samples 

Uncertainties in the e + e~ data sets are accounted for in several ways. For the e + e~ — > 7T°7 
data sets (351 [361 E7J and the e + e~ -» 777 data sets (381 EH [351 EH gOj, taking into account 
the magnitude of the systematics, we did not find motivated to split them up into their bin-to- 
bin correlated and uncorrected parts. We just add in quadrature the reported systematic and 
statistical errors. 

For the e + e~ — > samples, we dealt differently with the different data sets. For the 

relatively unprecise data sets [|4Tll42l we did as just explained for e + e~ — > (it /rj)^ data and 
simply added systematic and statistical errors in quadrature. Instead, for the more accurate data 
sets |fT9l 14311441 |45l , only the uncorrelated part of the systematic uncertainty has been added 
in quadrature to the statistical errors. On the other hand, the bin-to-bin correlated error has 
been treated as emphasized above for the r data sets. For reasons already emphasized in [QQ|, 
we have discarded the e + e~ — > data provided by 1146114711 . 

Finally, the various Novosibirsk e + e~ — > n + ix~ data sets recently collected [fT9ll20ll48ll2Tll 
carry a common bin-to-bin and sample-to-sample correlated uncertainty of 0.4 %. The older 
data from OLYA and CMD [49] also share a common correlated (scale) uncertainty [|50l of 

8 For instance, fitting the BELLE and CLEO spectrum lineshapes up to 1 . GeV does not reveal worse fit quality 
than when stopping the fit at 0.9 GeV. Higher mass vector meson influence in this region, if any, is thus found 
small enough to be absorbed by the other fit parameters. 
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~ 1%. In both cases, we have added the uncorrelated part of the systematics to the statistical 
errors in quadrature. 

Instead, in order to treat properly the correlated uncertainty, one should consider the data 
sets in Qj3|20l|48l|2T| as a single (merged) data set and use as x 2 an expression like Eq. (PT21) to 
introduce the common scale to be fitted. Here also, if A/0.4% < 1 -j- 2, one could remove this 
scale while performing the change V!y + (0.4%) 2 m i mj. One has performed the same 

way, mutatis mutandis, with the older OLYA and CMD data sets [49 1 . 

Because of the poor probability of the best fit to the KLOE data lfi~8l already commented 
upon in flU, the corresponding data sample is not included systematically in the set of e + e~ — > 
n + ix~ data samples considered ; however, its effects will be noted when relevant. In order to 
fit the parameters of the IB p, to, mixing scheme EMI], one still uses a subset of 9 radiative 
decay width data which are taken from the latest issue of the Review of Particle Properties E2l 
and are given explicitly in HI. 

3 Simultanous Fits To e + e~ and r Data 

3.1 Interplay Between 5m 2 , 5g And A 

Strictly speaking, the lineshape of the r spectrum is determined by the HLS parameters 
g' = g + Sg and the Higgs-Kibble mass m 2 + 5m 2 . This is clear from the expressions given 
in Section [2] above. The specific Isospin breaking parameters Sm 2 and Sg differentiate the p^ 
lineshape from those of the p° meson. However, these parameters also govern the absolute 
scale of the p^ spectrum compared to the p° one. Therefore, if an uncertainty on the absolute 
scale of a measured r spectrum calls for a fit parameter A rescaling the whole data spectrum, it 
is quite important to examine its interplay with Sg and Sm 2 . 



Data Set 


Sm 2 (GeV 2 ) 


Sg 


A 


X 2 / points] Exp 


ALEPH 


(3.37 ± 1.27) 10- 3 


(-0.56 ±0.12) 10- 1 


(-1.01 ±0.40)% 


27.16/38 


BELLE 


(-0.01 ±0.77) 10" 3 


(-0.12 ±0.10) 10- 1 


(-3.83 ±0.54)% 


32.46/20 


CLEO 


(-1.53 ± 1.07) 10~ 3 


(0.16 ±0.14) lO" 1 


(-5.51 ±0.74)% 


38.99/30 


ALEPH 


(4.04 ± 1.22) lO" 3 


(-0.69 ±0.11) 10- 1 





29.19/37 


BELLE 


(2.18 ±0.71) lO" 3 


(-0.51 ±0.08) 10- 1 





41.12/19 


CLEO 


(2.26 ±0.94) lO" 3 


(-0.54 ±0.11) 10- 1 





61.49/29 



Table 1 : Global fit results with each r data sample separately. Bolface numbers are actually 
not allowed to vary in the fits. Global fit probabilities are always above 90%. 



For the present exercise, we consider all data sets involving e + e~ — > it + tt~ data together 
with all data sets covering the e + e~ — > 7r°7 and 177 annihilation channels. These additional 
channels allow to remove the p°/co/(p — >■ (n° /rj)^ partial widths from the vector meson decay 
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mode subsample unavoidably used. In this Section, the ISR KLOE data sample for e + e — > 
n + n~ is removed from the collection of data sets to be fitted ; we also leave aside the e + e~ — > 
annihilation data which play a minor role in defining the vector meson mixing scheme. 

The e + e~ measurements with s < 1.05 GeV 2 are submitted to fit - in order to include the 
4> region - together with all r decay measurements with m n7r < 1 GeV. 

We have performed simultaneous fits of each of the A, B and C r data sets together with 
the e + e~ data referred to above. The results are shown in Table Q] and exhibit quite interesting 
features, depending on the particular r data set considered. 

The first three lines provide parameter values when 5m 2 , 5g and A are allowed to vary. For 
each r data sample, A is constrained by the relevant [\fit/f]Ex P ] 2 term in the global % 2 . In this 
case, one notes that : 

• The significance for a non-zero 5m 2 is at ~ 2.6 a for A and negligible for B or C , 

• The significance for a non-zero 5g is at the ~ 1 a level for for B or C but large for A 
(4.7a), 

• The values for some important correlation coefficients returned by the fit are large for 
each t data set : (8g, A) ~ (5g,5m 2 ) ~ —50 % and (A, 5m 2 ) ~ (25 50) %. These 
values reflect the interplay between 5g, 5m 2 and A in determining the absolute scales of 
the experimental spectra. 

• The significance for non-zero A's is data set dependent : 2.5 o\ for A, 7.1a\ for B 
and 7.40(Ta for C. Compared with the scale uncertainties induced by the errors on the 
respective B n7r , this corresponds to [\fit/rj\ A LEPH = 2.0 ± 0.8, [Xfu/v] belle = 2.5 ± 
0.35 and [Xfu/f]}cLEO = 3.2 ± 0.43. Taking into account the large correlations, between 
8g, 5m 2 and A, this looks to us acceptable. 

The corresponding fit residuals are shown superimposed in the upmost Figure [Q One 
clearly sees that the B and C residuals are well spread around zero. Those for A are slightly 
distorded around the p peak in a way opposite to B and C. The last data column in Table \T\ 
illustrates that each of A, B and C is well described by the global fit, simultaneously with e + e~ 
data (The so-called New Timelike data [US IM ED always yield x 2 /points ~ 118/127). 

The non-zero values for the A's could possibly mean that the reported measured values 
for B n7T are overestimated. However, the large correlations noted just above prevent such 
strong a conclusion, as the rescaling effect could well be absorbed by some different values 

for (Sg, 5m 2 ). 

It is, therefore, worth examining the case when A = is imposed. Of course, this turns 
out to prevent any rescaling of the experimental spectra and, thus, check wether the mass and 
coupling breaking we introduced could account alone for the measured reported normalizations 
of the various r spectra. The corresponding results are shown in the downmost three lines of 
Table [Hand call for important remarks : 

• The significance for a non-zero mass shift 5m 2 is always improved : 3.6 a for A, 3.1 a 
for B and 2.4 a for C, 
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• The coupling difference always becomes highly significant : the ratio of the central value 
to its uncertainty is 6.3 for A and B and 4.9 for CLEO. 

• However, fixing A = marginally degrades the BELLE description, more significantly 
the account of CLEO data, whereas the fit quality is unchanged - and quite good - for 
ALEPH data. 

The description of e + e~ data are, in this case, marginally degraded (The new timelike data, for 
instance yield x 2 — 125 instead of 118 before). 

This leads us to conclude that the additional Isospin Symmetry breaking mechanism, which 
summarizes into mass and width differences for the p's, allows to account for the original 
absolute normalization of the ALEPH spectrum with a very good probability. One can thus 
conclude that the dipion ALEPH spectrum - or, equivalently, the ALEPH pion form factor - 
is in in full accord with VMD predictions, provided one suitably accounts for mass and width 
shifts in the p ± information compared to p°. Numerically, this turns out to plug into the ALEPH 
spectrum parametrization : 

f I >^> ^ ™l = [-2.50 ± 0.40] % 

I r p° 9 (13) 

[ m 2 p± - mj = 6m 2 = (4.04 ± 1.22) 10~ 3 GeV 2 

which approximately^] means m p ± — m p o ~ 2.59 ± 0.78 MeV and T p ± — T p o ~ —3.7 ± 0.6 
MeV. If 5m 2 , is numerically in the usual ballpark (|22|. Sg seems slightly larger than could have 
been inferred from commonly reported estimates [|22l l23l for T p ± — T p o . 

The results corresponding to fits with A = are shown in the downmost Figure [Q Com- 
pared with the upmost one, this Figure shows interesting features : the ALEPH residual dis- 
tribution is essentially unchanged, but the residuals for BELLE and CLEO are shifted towards 
positive values and start resembling the ALEPH residual distribution. This indicates graphi- 
cally that some rescaling is still needed in order to get a satisfactory description of the B and C 
data sets. Whether this residual rescaling could be absorbed into a more sophisticated Isospin 
breaking distortion mechanism cannot be discarded. 

As a summary, one can assert that, with appropriately chosen 5m 2 and Sg, ALEPH data 
can avoid any significant rescaling within the model we developped. This attractive property 
is unfortunately not completely shared by BELLE and CLEO in isolation. This is reflected 
by the fit information given in the last data column of Table \T\ and by the downmost Figure 
[0 Nevertheless, with some rescaling, the BELLE and CLEO data sets can be satisfactorily 
understood. 

Moreover, it is clear from the above analysis that there is a large interplay between param- 
eters defining IB shape distortions and the absolute scale of the experimental spectra. Stated 
otherwise, some part of a fitted rescaling could be due to the difficulty to model the actual shape 
distortions and, conversely, some part of a real scale factor could well be absorbed by a fitted 
distortion. Therefore, one cannot claim, in view of a large fitted rescaling, that the reported 
B n7r should be rescaled by as much. 

9 Mass and width of broad objects like the p are, conceptually, definition dependent and, additionally, we have 
not accounted for the large negative correlation term (5g, 5m 2 ). Moreover, these numerical values somewhat vary, 
depending on the exact e + e~ data set content submitted to the global fit. 
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3.2 Fitting The r Lineshapes 



In view of the different behaviors of the A, B and C data sets, it is worth examining sep- 
arately the r spectrum lineshapes besides the full spectra. As explained in Subsection I2.2.2L 
this turns out to fit the r data samples without including the scale constraining terms in their 
contributions to the total x 2 ■ 

For this purpose, one has first performed a global fit leaving free Sm 2 , Sg and the A's 
unconstrained. One thus gets Sg = (0.15 ± 0.08) 10" 1 and 5m 2 = (-1.47 ± 0.67) 10~ 3 GeV 2 
witr@ Xaleph = 17-86, xIelle = 29.37 and xIleo = 29-71 associated with a X 2 = H8.31 
for the 127 data points of the so-called New Timelike data |fi~9l l20l l48l 1211 . Instead, when 
imposing 5g = 5m 2 = 0, one gets Xaleph = 19-22, Xbelle = 27.61 and Xcleo = 3L83 
while the New Timelike data get x 2 = 120.23. 

One may compare the x 2 values yielded for each of the r data samples when fitting together 
their normalized spectra with those already reported^ in Table [TJ 

The results given just above enforce the conclusions we reached in Subsection 13. II In the 
global fit procedure, it looks difficult to distinguish effects due to IB shape distortions in the 
p distribution from genuine rescaling effects. This is especially striking when considering the 
ALEPH data set : One gets a quite acceptable description of the invariant mass distribution by 
either assuming A = and (Sg, 8m 2 ) free or by letting A free and imposing Sg = 5m 2 = 0. 
So, the reported scale discrepancy can be absorbed ; however, the sharing of the effect between 
true (physical) lineshape distortions and an (experimental) bias in the accepted value for B n7T is 
beyond the scope of fit results. This is an important issue for using r data in order to estimate 
hadronic contributions to g — 2 ; this will be commented on at the appropriate place. 

Therefore, at least when fitting the r lineshapes, it is not justified to keep non-vanishing 
Sm 2 and Sg. In this case, the main difference between the p + and the (physical) p° is essentially 
carried flU CO by the 7 — p° and the W — p ± transition amplitudes f< and f T p : 

| = l + ^(s) + ^/3(s) , U; = ag 2 fl) (14) 

These amplitudes differ through terms explicitly depending on the Isospin Symmetry break- 
ing functions which account for the vector meson mixing [@] [fl . The s-dependent terms in 
Eq. (fl"4l) account for the isospin zero component of the p° meson. Keeping non-vanishing 
5g and Sm 2 within our modelling, the right-hand side of Eq. (fT4"l) gets a leading additional 

(8m 2 jm 2 — Sg/g) term. 



From now on, we choose to reintroduce the e + e~ — > data from PTTI 1421 and 

lfT9l l43l l44l l45l . as discussed in Subsection 12.2.31 above. We also let the scale factors vary 
for both the e + e" and r data sets^l. Some numerical results obtained with various data set 
configurations are reported in Table [2l One notes that all scale factors introduced for the 
e + e~ — > n + ix~ and e + e~ — > cross sections are in good correspondence with the 

10 The number of data points in the r data samples submitted to fit are always 37 (A), 19 (B) and 29 (C). 
n Let us remind that the fit results reported in this Table are, instead, obtained from (global) fits using each of 
the r data sets in isolation. 

12 keeping, of course, the term constraining the scale variation for all the relevant data sets. 
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Expected 


No r E 


)ata 


t Data Set ( 


Configurations 




r.m.s. 


e+e" NSK+KLOE 


e+e" NSK 


NSK + A sh B sh C sh 


NSK + ABC 


Scale New Timelike 


0.4% 


(-0.9 ±0.4)% 


(-0.6 ±0.4)% 


(-0.6 ±0.4)% 


(±0.4 ±0.3) % 


Scale Old Timelike 


1.0% 


(+1.4 ±0.7)% 


(±1.6 ±0.7)% 


(±1.6 ±0.7)% 


(±2.4 ± 0.7) % 


Scale CMD-2 El 


1.3% 


(-0.4 ± 1.2)% 


(-0.5 ± 1.2)% 


(-0.4 ± 1.2)% 


(-0.7 ± 1.2) % 


Scale CMD-2 flU 


2.5% 


(-2.6 ±1.9)% 


(-1.8 ± 1.9)% 


(-1.8 ± 1.9)% 


(-1.4 ± 1.9) % 


Scale CMD-2 flU 


4.6% 


(-4.7 ±3.4)% 


(-4.2 ± 3.4)% 


(-4.1 ±3.5)% 


(-3.6 ± 3.4) % 


Scale CMD-2 flU 


1.9% 


(-2.6 ± 1.6)% 


(-2.2 ± 1.6)% 


(-2.0 ± 1.6)% 


(-2.1 ± 1.6) % 


9 


— 


5.569 ±0.003 


5.587 ±0.008 


5.565 ±0.006 


5.532 ±0.007 


5g 


— 


— 


— 





[-0.30 ±0.07] lO" 1 


5m 2 (GeV 2 ) 


— 


— 


— 





[1.09 ±0.6] 10~ 3 


X 


— 


0.917 ±0.012 


0.908 ±0.013 


0.905 ±0.013 


0.902 ±0.013 


ZA 


— 


1.501 ±0.010 


1.476 ± 0.009 


1.468 ±0.009 


1.454 ± 0.009 


a 


— 


2.372 ±0.002 


2.357 ±0.005 


2.361 ±0.004 


2.382 ± 0.004 


c 3 




0.927 ±0.006 


0.936 ±0.006 


0.944 ±0.006 


0.952 ±0.006 


Cl - c 2 




1.194 ±0.031 


1.196 ±0.032 


1.228 ±0.032 


1.237 ±0.032 


X 2 /dof 
Probability 




647.43/653 
55.4% 


521.15/593 
98.5% 


602.36/675 
97.9% 


645.22/676 
79.7% 



Table 2: Results in various fit configurations. e + e~ NSK data stand for all annihilation pro- 
cesses discussed in the text. Inclusion of ALEPH, BELLE and CLEO data are referred to as 
A, B and C spectra respectively, A sh , B sh and C sh denote the corresponding normalized spec- 
tra. The upper part displays the rescaling factors for n + ix~ (first two lines) and data 
samples. The lower part refers to breaking parameters discussed in Ref. OX Numbers written 
boldface are parameter values not allowed to vary inside fits. 

expectations reminded in the first data column. This result is understood as confirming once 
more [1J the claims of the corresponding experiments concerning their correlated systematic 
uncertainties. In the same Table, the data column A s/l B s/l C s/l reports the results obtained while 
removing the scale constraining terms (the second one in Eq. (PT21) for each r data sample). The 



13 



numerical values of the rescaling parameters for the r spectra are given below. 

Table[2]justifies the removal of the e + e~ rescaling factors^ from the set of parameters to be 
fitted. One also observes that the physics parameters have only (small) reasonable fluctuations. 
Finally, Table [2] clearly shows that all Novosibirsk (NSK) data (4 different cross sections) are 
quite consistent with the freely rescaled r decay data as the fit yields a probability above the 
90% level. Moreover, it is also interesting to note that constraining the rescaling of the r 
spectra by their reported scale uncertainty (i.e., keeping the scale fixing terms as in Eq. (fT2l) ) 
still provides quite comfortable probabilities (above the 80% level), as can be seen from the 
last data column in Tabled 

This last data column exhibits also interesting features : In a simultaneous fit to all r data 
samples together with e + e~ data and vector meson decay modes, the significance for a non- 
vanishing 5g is ~ 4er whereas that for 5m 2 is only ~ 1.8er. Comparing with the separate r data 
sample fits reported in Table Q] indicates that significantly non-vanishing Sm 2 and 5g are driven 
by ALEPH data only. 

We have explored the effect of including the KLOE data set [[TBI beside the whole set of 
Novosibirsk e + e~ data and the three r decay data sets. As expected from its intrinsically large 
X 2 ED (123 for 60 measurements in the present case), the global fit probability drops to ~ 8%. 
However, the information displayed in Table |2] is not substantially modified. In particular, the 
scale factors affecting the e + e~ Novosibirsk data remain very close to the expectations reported 
in the first data column. 

From now on, all scale factors affecting the e + e~ data - except for KLOE, when relevant 
- are set to zero and the scale uncertainties are transfered into the error covariance matrix as 
explained above and emphasized in Section 6 of [HI. 

3.3 Final Global Fits To e+e and r Data 

In this Subsection, we only refer to fits performed without the KLOE data set. Some results 
are displayed in Table [3] and will be discussed now ; other pieces of information will be given 
later on, especially those with each r data sample fitted in isolation with the whole set of e + e~ 
data. 

Two kind of results are shown in Table |3l the former with fitting only the r spectrum 
lineshapesL^ (see Subsection 12.2.21 above), the later with the full r spectrum information as 
expressed by Eq. (fTTT) . The first data column in Table |3] displays the results obtained with only 
the whole set of e + e~ data d). 

The first two data columns in Table[3]clearly illustrate that the fit quality remains optimurrfjfl 
(above the 90% level) when including only the r spectrum lineshapes inside the fit procedure 
(data column flagged with A sh B sh C sh ). This does not hide some badly described data set as 
reflected by the various partial x 2 displayed. Figure © displays the normalized spectra together 
with the fit function ; the residual distributions shown in Figure © confirm the goodness of the 
fit. 

13 Then, the correlated scale uncertainties only appear in the full error covariance matrix as reminded in Sub- 
section |223] and in Section 6 of QJ. 

14 In this case, the t data are flagged as A , B sh and C sh for respectively the ALEPH, BELLE and CLEO data 
sets. When using, for each t data set, the full Eq. il2i . these data sets are simply flagged A, B and C. 

15 The origin of such favorable probabilities is discussed at the end of this Subsection. 
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Comparing the individual x 2 m the second data column with those in the first one for the 
different e + e~ data subsets (and the decay width data set), one clearly observes negligible 
modifications of their fit quality. Moreover, the x 2 obtained for the individual A sh , B sh , C sh 
spectrum lineshapes are also quite reasonable. 

The numbers shown boldface and within parentheses in the first data column are also quite 
interesting. Fixing the arbitrary scales to their fit values (given in second data column) for each 
r data set, one can compute the x 2 distance of each of the A sh , B sh , C sh spectrum lineshapes 
to the fit solution to only the e + e data. Comparing these numbers to their homologues in the 
second column (obtained, instead, with the r spectrum lineshapes fitted) clearlyproves that the 
e + e~ data allow to predict the r spectrum lineshape with very good accuracy^. One should 
note that, if unconstrained, all r spectrum normalizations prefer consistently a rescaling by 
~ —5%. As discussed in Subsection l3.1l above. the exact meaning of such rescalings should be 
considered with great care. 

In Figure |4] we plot the fit residuals normalized to the fit function value, namely : 

(l + X)mi-fi . 
Xi = , i = l ) N Exp (15) 

Ji 

neglecting the effects of correlations on the plotted errors. This plot clearly shows that, up to 
~ 1 GeV, the residuals can be considered as flat, structureless, distributions ; Figured is the 
exact analog of Figure 12 in ifTTl and exhibits the similar character of BELLE and CLEO data 
in our fitted s range. 

The last data column in Table [3] displays the fit information while fitting the B^/N dN/ds 
distributions of the various r experiments together with the whole set of e + e~ data, taking 
into account the constraint imposed on the A scale factors by the uncertainty on the measured 
branching ratio values B^. 

The full fit thus obtained is also quite satisfactory as reflected by the fit probability (80%). 
It shows, nevertheless, that a significant rescaling of the experimental data is requested. For the 
r data sets, the partial x 2 information is displayed as x/y where x is the part of the x 2 coming 
from the data points and y — x is the contribution of the (A/77) 2 term ( see Eq. (fT2l)). y/y — x 
tells how far from the B^ central value - in units of r/Exp - , the spectrum normalization is pre- 
ferred. Assuming no external cause to the rescaling, the numerical values of the corresponding 
coefficients look close to the expected B nn value, taking into account the reported accuracy of 
the various measurements. Compared to these, the fit values are found at resp. 2.9 (Jaleph, 
1.4 ubelle and 1.5 ctcleo towards lower central values for B nn . 

The results are represented in Figures ©, © and ©. In Figure © one clearly sees that the 
residuals for CLEO and BELLE are consistent with those in Figure ©, while those for ALEPH 
exhibit a structure around the p peak location. We have widely discussed in Sub section [3J] the 
behavior of ALEPH, BELLE and CLEO residuals under very close conditions^- Figure © 
gives another view : Comparing Figure © to Figure © shows that the residual distributions 
remain flat up to the ~ <p mass, but slightly shifted towards positive values. 

16 We already reached this conclusion |4| with the CLEO data sample lineshape. 

17 The differences between Figures|6]and[T]is essentially due to the fact that we now perform a simultanous fit of 
ALEPH, BELLE and CLEO data sets instead of examining each them in isolation. In these conditions the values 
found for 5g and 8m 2 are dominated by BELLE and CLEO which prefer a smaller value for Sg than ALEPH. 
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Data Set 


No r Data 


t Data Set ( 


Configurations 


(fj data points) 


NSK(e+e-) 


NSK +A sh B sh C sh 


NSK + ABC 


Decays (9) 


14.84 


15.12 


16.20 


New Timelike (127) 


118.88 


120.24 


126.47 


Old Timelike (82) 


50.65 


51.08 


60.45 


7T°7 (86) 


66.03 


65.84 


66.07 


?77(182) 


135.26 


135.54 


135.78 


TT + TT-TT (126) 


139.45 


139.42 


139.44 




— 





[-0.30 ±0.07] 10- 1 


5m 2 (GeV 2 ) 


— 





[1.06 ±0.62] 10~ 3 


ALEPH (37+1) 


(23.52) 


19.23 


28.30/36.51 


ALEPH Scale (%) 


— 


-5.45 ± 0.62 


-1.46 ±0.38 


CLEO (29+1) 


(35.57) 


31.83 


37.22/39.46 


CLEO Scale (%) 




-5.48 ± 1.01 


-2.60 ±0.54 


BELLE (19+1) 


(27.16) 


27.61 


26.43/28.29 


BELLE Scale (%) 




-4.80 ±0.71 


-2.09 ±0.46 


X 2 /dof 
Probability 


525.10/597 
98.4% 


605.90/679 
97.9% 


648.68/680 
80.1% 



Table 3: Individual x 2 of the data samples under various fit configurations. The e + e~ data sets 
are referred to in the text. Inclusion of ALEPH, BELLE and CLEO data sets are referred to in 
column flags as A, B and C, respectively. A sh , B sh , C sh correspond to fitting the lineshape of 
the normalized invariant mass spectra only. Numbers indicated boldface are the x 2 distances 
of A sh , W h , C sh to the fit solution of only the e + e~ data. Parameter values written boldface are 
not allowed to vary in fits. 

The error covariance matrix of the fit parameters still exhibits some large correlations as 
(g, Sg) ~ 40% and (g, 5m 2 ) ~ (Sg, 8m 2 ) ~ —40%. One also gets correlations at the +20% 
level between Xaleph, ^belle and Xcleo, while all others do not exceed a few percent level. 
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The sign of the correlations for (g, 8m 2 ) and (8g, 5m 2 ) prevents a possible relation of the 
form m 2 + 5m 2 = a[g + 8g] 2 f 2 for the p^ mass squared. Relying on the global fit of all r data 
samples, one gets : 



significantly smaller than for ALEPH only (see Eq. ([T3T)). In this case, the mass difference 
between the two p mesons is m p ± —m p o ~ 0.68±0.39 MeV, whereas T p ± —T p o ~ — 1.63±0.39 
MeV only. 

Summarizing the results gathered in Tabled one can conclude that e + e~ data and r data 
do not exhibit inconsistencies. Indeed, the fit of the spectra reported in the third data column 
is quite satisfactory, and this takes into account all reported experimental information. This 
agreement at the the form factor level, does not a priori mean that similar g — 2 values will be 
obtained when using/removing r data. 

Taking into account the statistical significance of the rescaling factors as reflected by their 
uncertainties (see Table |3]), we neglect none of them in the rest of this analysis. One may 
stress once again that these rescaling factors may well only reflect an incomplete account of 
the isospin breaking effects affecting the r spectra lineshape which have to be "subtracted" in 
order to compute estimates of the nn loop contribution to g — 2. 

Finally, one may also wonder to reach as frequently such favorable probabilities. As can be 
seen in Tabled this reflects the low contributions to the x 2 yielded by the old irn data ll49ll5TI 
( x 2 /n ~ 0.6) and by the (recent) e+e" -»■ (tt%)7 data GIC|33imi2ilE71ED (X 2 /" ^ 0.75). 
Instead, the newly collected nix data sets [fP9l |20l |2TT| yield x 2 / n — 0.9 and the 3-pion data 
[[191 |43l Sll |45l ST] g2]| X 2 /^ — 11, quite comparable to CLEO data x 2 /n ^ 1.3. Discarding 
the contribution of the scale penalty terms, BELLE data yield x 2 / n — 1-4 and ALEPH ~ 0.9. 
Therefore, these high probabilities reveal certainly some overestimation of systematic errors in 
definite e + e~ data sets. When data sets for e + e~ — > (w° /rf)^ with better estimated errors will 
become available, this question will be naturally solved ; awaiting this, one may conclude that 
the most precise pieces of information provided presently by the full e + e~ — > (n° /rj) f y cross 
sections are essentially the partial widths for p°/co/(p — > (ir° /rj)^. 

On the other hand, this explains why one should carefully control that some inconsistency 
in the fit is not hidden by these favorable probabilities. A good criterium is certainly provided 
by the x 2 l n value associated with each data set and by its behavior while modifying the data 
set collection submitted to fit. 

4 Consequences for the Hadronic Contribution to g — 2 

An important outcome of our model and of our treatment of the data sets - i.e. our fitting 
procedure - is related with the estimate of hadronic contributions to g — 2 up to 1 GeV. Mixing 
different processes correlated by the same underlying physics cannot be successful without 
some clear understanding of the errors in each data set, which should be properly implemented 
inside the fitting code. Our procedure relies on the whole available information on uncertainties 




(1.06 ± 0.62) 10~ 3 GeV 2 



[-1.10 ±0.26] % 



(16) 
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(magnitude and type). This allows us to draw conclusions directly from a global fit to several 
data sets rather than combining physics information derived from separate fits to the individual 
data samples. 

However, this does not guarantee that combining various experiments will result in im- 
proved uncertainties. Indeed, the way the systematic errors, especially scale uncertainties, 
combine in the fit procedure cannot be straightforwardly guessed. One may, nevertheless, ex- 
pect the results to be more reliable ; indeed, assuming systematics are randomly distributed, 
the net result of mixing different experiments and/or data sets should be to neutralize them to a 
large extents 

The lowest order contribution of a given annihilation process e + e~ — > H to the muon 
anomalous magnetic moment a M = (g — 2)/2 is given by : 



" ,,{H) = fr ds k ^ s) ^ a7) 



where a(s) is the Born cross section of the annihilation process e + e~ — >■ H, sh its threshold 
squared mass and s cut is an assumed end point of the non-perturbative region. K(s) is a known 
kernel [fT3l given by the integral : 

K(s)= fdx 2 , (18) 

Jo x + (1 — x)s/m^ 

m M being the muon mass. For s > &m 2 u , this writes : 

In (1 + x) — x + 



K(s) = X 4(2-x>) + ( 1 + x2 ^ + x2) 



2 y ' x 2 



x 2 

T 



1 + X 2 

H x mx 



1-/3 / Am 2 
with : x = and 8 = \ 1 

1+8 H V s 

and, for < s < Arn 2 ^, it becomes lfT3l (r = s/m 2 ) : 



(19) 



K(s) = - - r + -r(r - 2) lnr - (1 - 2r + — )J—^— arctan J- — - (20) 
2 2 2 y 4 — r V r 

This expression has to be used in order to integrate the cross section for e + e~ — > 7r°7 below 
the two-muon threshold. 

Our global fit provides the theoretical Born cross sections with their parameter values, 
errors and their (full) covariance matrix. As illustrated above and in [1J, the results obtained 
while fitting altogether several e + e~ cross sections and r spectra are satisfactory. Therefore, 
we consider that using our cross sections within a Monte Carlo, which fully takes into account 
the parameters, their errors and correlations, should provide a fairly well motivated value for 
each accessible a^(H) in our fitting range {i.e. from each threshold up to ~1 GeV/c) and for 
its uncertainty which appropriately merges statistical and systematic errors. 



18 Of course, one cannot completely exclude that systematics could "pile up" coherently ; however, if the 
number of different data sets is large enough, such a possibility looks rather unlikely. 
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In order to allow for a motivated comparison between our results and experimental esti- 
mates for a A1 (7r + 7r _ ), one should also include the effects of Final State Radiation (FSR) into 
our estimates of the ixix contribution to the muon anomalous moment. Indeed, even if not al- 
ways manifest, the shift produced by FSR corrections is included in the reported |fT9l |20l I2T1 
experimental values for a^(T{ + n~). The FSR corrections are accounted for by multiplying the 
7T7T Born cross section in Eq. (fT71 ) by |fT3l : 

1 + Vfsr(s) = (l + v(s)^r ~ ^) ^ (l - exp {~^}\ (2D 

where the Schwinger function ll52ll r](s) can be found corrected for a missprint in ll53l together 
with a simplified expression (also derived by Schwinger) valid at the 1.5% level. The uncer- 
tainties affecting the FSR effect estimates are not known ; nevertheless, they are expected small 
in our range of interest [fT3l ; we, therefore, neglect their contribution to the errors we report 
on a M (7T7r). Finally, FSR effects on contributions other than txtx to a M are also known to be 
negligible lfT3l and are thus neglected. 

4.1 Global Fit Results And e+e Data 

It is quite important to compare the outcome of our model and fit procedure with experi- 
mental data. This should allow to check for possible methodological biases and to substantiate 
how the merging of statistical and systematic errors operates. 

CMD-2 QH [201 and SND ED, for instance, have published a^ix+TX') obtained by in- 
tegrating numerically their measured e + e~ — > cross section over the interval y/s = 
0.630 -T- 0.958 GeV, using the trapezoidal method. We have run our code using separately each 
of the data sets given in resp. lfT9l . EDI and f2\\ together with the relevant set of decay par- 
tial widths of vector mesons needed in order to determine numerically the SU(3)/U(3)/SU(2) 
breaking parametrization. 

The results are given in the upper part of Table |4] and look in good correspondence with ex- 
pectations. The highly favorable probability for the SND spectrum 11211 might be attributable 
to its larger systematics compared to CMD-2. The middle line in this Table shows the re- 
sult obtained while fitting simultaneously the three Novosibirsk spectra |fT9l l20l |2TI together. 
The global fit of these data sets, where systematics are certainly under good control, provides 
a strong reduction of the global uncertainty - merging all reported systematic and statistical 
errors. The improvement of the uncertainty derived from the global fit solution compares fa- 
vorably with the average value proposed by |fl4l using a spline method. Therefore, nothing 
obviously abnormal is recognized in this comparison. The following line shows that KLOE 
data [54] allows to reduce a little bit more the uncertainty ; both the central value and the un- 
certainty compare well with the average proposed by lfl4l . However, one should note the very 
low probability of the associated global fit. 

The next three lines in Table |4] are also quite interesting. The first of these reports the result 
obtained when fitting using all available data sets - except for KLOE |[T8l ; this turns out 
to include into the fitted data the so-called "Old Timelike Data" [|49l IBTI (82 measurements) 
besides the "New Timelike Data" (127 measurements). The central value for a^nn) is slightly 
shifted downwards by about la with a negligible gain in precision. One also observes that the 
fit probability makes a large jump (~ 50% 90%), reflecting the larger uncertainties in 
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Data Set 


Experimental Result 


Fit Solution 


Statistical Information 








X 2 /dof 


Probability 


CMD-2 (1995) 


362.1 ± (2.4) stot ± (2.2)^ 


362.57 ±2.64 


42.46/44 


53.8% 


CMD-2 (1998) 


361.5 ± (1.7) stat ± (2.9)^ 


362.36 ± 2.14 


38.05/40 


55.9% 


SND (1998) 


361.0 ± (1.2) stat ± {A.7) syst 


361.09 ± 2.04 


27.10/46 


98.8% 


NSK New Timelike 7r + 7r~ 


360.0 ± 3.02 exp *** 


361.39 ± 1.72 


124.61/128 


56.8% 


\ i p i/ \t rp" 1 • 1 _1_ — „ 1 S 1 /^\ 1 — 1 

NSK New Timelike 7r + 7r + KLOE 


or'op" i nil 

358.5 ± 2.41 exp 


360.25 ± 1.48 


OCT r r)Q /inn 

ZOO.ZZ/ loo 


f\ 1 07 


INolV^Tr 7T J 




qccq t^n -i- i fin 
ooy.ou ± i.ou 


180.09/210 


93.3% 


NSK(7r+7r-) + (7T°/r/)7 




359.42 ± 1.52 


373.59/468 


99.96% 


NSK (tt+tt-) + (7r°/r/)7 + (tt+tt-tt ) 




359.31 ± 1.62 


525.10/597 


98.4% 



Table 4: Contributions to 10 10 a M (7r7r) from the invariant mass region 0.630 — 0.958 GeV/c. The 
data flagged by *** are combined values proposed in [14] for the New Timelike data altogether, 
with or without KLOE data [|54l . The first three lines in this Table refer to a^nn) experimental 
values given in resp. JTH, fl20l and BUI. "NSK (7r + 7r~)" indicate that all annihilation data to 
7r + 7r~ are considered. 



the data from ll49l IBTI compared to those from lfi~9l l20l |21~T| . This once more shows that the 
systematics in the data sets from [491I5TI have been conservatively estimed. 

The last 2 lines in Table |4] exhibit the same trend ; indeed, including in the fit procedure the 
86 measurements of e + e" — > 7r°7, the 182 measurements of e + e~ — > 777 and the 126 points 
of e + e~ — > from |fT9l |43l 1441 |45l I4T1 1421 does not improve, strictly speaking, the 

uncertainty ; however, it is satisfactory to check that the central value only fluctuates inside 
quite acceptable limits. This teaches us some remarkable facts : 

• Improving systematics in all processes having the same underlying physics as ir + n~ 
might be useful to allow a better determination of their own contributions to a M , but also 
of a^irn) itself. Indeed, most physics parameters in the quoted processes are the same 
as in e + e~ — > n + ix~ . Conversely, higher quality n + n~ data should improve estimating 
the contributions of the other related annihilation processes to a M . 

• If errors are reasonably well understood and appropriately dealt with, adding poor qual- 
ity data sets into the fitted sample does not degrade the result ; this allows, however, to 
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confirm the stability of the central values, which is certainly a valuable piece of informa- 
tion. 

However, the most important remark is certainly that using the radiative partial width de- 
cays of light mesons - and/or the e + e~ — > (it /if)^ cross sections - together with e + e~ — » 
n + n~ allows a quite significant improvement of the accuracy on a^nn) compared to the di- 
rect numerical integration of the experimental spectra. This is, indeed, the main advantage of 
having a constraining global model allowing for an overconstrained global fit. 

4.2 Effects Of Including r Spectra And KLOE Data 

One may question the effects produced by introducing the r decay data inside our collection 
of fitted data sets. These are expected to improve the model parameter values (and errors), if 
their systematics are indeed reasonably well controlled. 

Our strategy will be to consider all e + e~ data sets used just above and look at the effect of 
including the A [|T5l . B ifTTl and C lfT6ll data samples in isolation or combined. 

As above, when using the fit results obtained by constraining the rescaling factors - i.e. 
keeping the (A/77) 2 terms in the x 2 expressions for r data samples - the r samples will be 
denoted A, B and C ; when concentrating over r lineshapes, these will be denoted A sh , B sh and 
C sh . The case when fitting the r spectra by allowing non-vanishing 5g and Sm 2 and imposing 
A = will still be referred to as Ad m ,d g ^dm,d g or Cd m ,dg in our Tables and/or Figures. 

4.2.1 Comparison With Standard r Based Estimates for a M 

When letting free 5g and Sm 2 while imposing A = 0, our approach is comparable to those 
underlying the so-called r based estimates lfl4l of the hadronic contribution to the muon g — 2. 
Indeed, the whole set of e + e~ data fixes in a data driven mode all identified isospin breaking 
corrections : p — cu and p — <p meson mixing, p meson mass and width differences, 1=0 part of 
the p revealed by its coupling to photon0 (see Eq. ([141) ) and FSR corrections. The pion mass 
difference is plugged in directly. It is thus interesting to compare our results in this case with 
the existing estimates and thus check the effects of a global fit. The hadronic contributions 
to the muon g — 2 derived from r data have been updated recently in [14] ; the numerical 
contributions to from the the reference region y/s £ [0.63, 0.958] are not published but have 
been kindly communicated to us 11551 . 

Table[5]displays the experimental data derived from each existing r data set and their com- 
bination Ifl4ll55l in the first data column. One may note that the proposed experimental average 
is larger than each of the individual estimates ; actually, in order to perform the average, each 
individual estimate has been rescaled to the world average value for B nn ll55l . The total (exper- 
imental) errors displayed are computed by summing up the various components in quadrature. 
This is provided in order to allow for an easy comparison with our fit result. The data columns 
X^/dof and xL/dof display the contribution to the total x 2 provided by resp. the r and the 
so-called New Timelike e + e~ annihilation data |[T9l 12011271 . which serve as quality tags. 

19 This contribution is currently not considered as such lfT4ll . It should be partly absorbed in p meson width 
corrections. 
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Data Set 


a M (7T7r) 




Statistical Information 




Experimental Result 


Fit Solution 


X 2 T /dof 


Xee/ dof 


Probability 


A dm , dg 031 


364.02 ± (2.19) exp ± (1.97) flr ± (1.51) JB 
(364.02 ± 3.31 tot ) 


362.83 ± 1.46 


45.52/37 


122.12/127 


95.3% 


Bdm,dg El 


366.44 ± (1.02) exp ± (5.70) Br ± (1.51) /B 
(366.44 ±5.98 tot ) 


364.85 ± 1.32 


41.95/19 


128.58/127 


76.5% 


Cdm,dg EH 


366.62 ± (A.17) exp ± (6.37) Br ± (1.51) Ifl 
(366.62 ± 8.05 tot ) 


364.23 ± 1.82 


63.01/29 


125.63/127 


70.0% 


OPAL (2D 


354.40 ± (4.67) e:cp ± (4.78) Br ± (1.51) /B 
(354.40 ± 6.85«) 










ALL r Sets 


367.46 ± (1.31) e2p ± (1.59) Br ± (1.51) /B 
(367.46 ± 2.55**) 


367.12 ± 1.30 


102.06/85 


143.32/127 


58.2% 



Table 5: Contributions to 10 10 a M (7r7r) from the invariant mass region 0.630 — 0.958 GeV/c. The 
experimental values If55l are derived using the method from Ifl4l . The experimental average 
includes OPAL data, our fit result does not. The meaning of Ad m ,d g , Rdm,dg> C-dm,dg is explained 
in the text. 



One should note the nice correspondence between our central values and the corresponding 
experimental estimates. The improvement of the total errors provided by the global fit method 
is also worth mentioning. The reduction of the uncertainties provided when using r data looks 
even more important that when using e + e~ data alone. Of course, the errors provided there, as 
anywhere in this paper, are the minuit errors returned by the fits. 

The result for (Ad m ,dg ^dm,d g Cd m ,dg) - last line in Table |5]- is also quite remarkable. It 
clearly illustrates that our global fit does not lead to a standard averaging, but takes into account 
the relative mismatch of the A, B and C lineshape distortions noted in Subsection 13.11 In this 
procedure, the fit average is significantly pushed upwards, in accord with the experimental 
estimateo 

The data columns providing the x 2 information are also quite important. As already noted 
in Subsection 13 . 1 L Xaleph i s reasonably good simultaneously with the New Timelike e + e~ 
data. Comparing Xaleph h ere (45.52) with its homologue in Tabled] - fourth line therein 
- (29.17), reveals that the ALEPH data sample meets some difficulty in accomodating the 
data for e + e _ — > However, the description remains quite reasonable and one may 

conclude that there is no real mismatch, within our approach, between ALEPH data and the 
VMD expectations. The value just derived from ALEPH data compares well with those 

20 This should be slightly larger, if removing the OPAL data set J3T|. 
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derived from e + e~ data only (see Table HI). 

The values for a M derived for BELLE and CLEO data in isolation, even if slightly larger 
than expected from VMD inference, are not in strong disagreement with these. However, as 
can be seen from either of Tableland Table[T|(see the A = entries therein), the fit quality, as 
reflected by Xbelle an ^ Xcleo l°°ks significantly poorer (both yield x 2 '/npoints ~ 2) than 
ALEPH (xV^points ~ 1.1 1.2). Whether this behavior reveals specific systematics is an 
open issue. 

When performing the global fit with all r data samples (last line in Table [5]), the weight 
of r data happens to become dominant. In this case, the fit of the r data is nearly unchanged 
(Xaleph = 26 - 64 ' X 2 belle = 45 - 89 > X 2 cleo = 29 - 52 ) I however, even if apparently reason- 
able, the new timelike e + e~ data yield a x\ e increased by 20 units, which is significantly far 
enough from optimum that we consider the corresponding a M value with great care. 

The analysis in Subsection 13.11 has shown that rescaling factors and lineshape distortions 
are sharply related. However, in order to derive r based estimates for a M , one should remove 
all distortions produced by isospin breaking effects specific of the meson relative to p°. 
Therefore, we do not consider reliable the estimates provided in Tabled except possibly those 
in the A dm ^ dg entry, as the distortions in the B and C samples are to be better understood. 

4.2.2 Additional r Based Estimates of 

Our analysis in Subsection 13. 1 1 provided a serious hint that rescaling the r spectra can be 
an effective way to account for shape distortions produced by IB effects which differentiate 
the p^ 1 and p° mesons. We just provided results assuming no scale correction to the r data 
samples {Xaleph = Xbelle = Xcleo = 0). In order to provide a complete picture, it is 
worth considering two more sets of fitting conditions. In this Subsection, we will present the 
results derived by : i/ fitting the r spectrum lineshapes by the method emphasized in Subsection 
12.2.21 nl fitting the r spectra by allowing non-vanishing rescalings, however, constrained by 
the relevant (X/r) Exp ) 2 term for each of the r data set. 

The upmost part of Table [6] collects our results for a M obtained using the r data samples, 
except for the first line which is nothing but the final result in Table HI The second line shows 
that including the r lineshapes into the fit (A sh , B sh , C sh ) gives results perfectly consistent with 
using only e + e~ data and a slightly improved uncertainty. As the fit quality of each A sh , B sh 
and C sh on the one hand, and each of the e + e~ data are simultaneously optimum, one may 
consider this estimate safe. 

Compared with pure e + e~ estimates, the next 3 lines show the effects of including each of 
A, B and C in isolation^. The value found for ALEPH agrees with that in the Ad m ,dg entry 
in Table [5J Interestingly, the values found using either B or C are significantly lower than their 
homologues in Table[5J but in nice correspondence with VMD expectations (see Table[4]). This 
may indicate that some rescaling is needed for these, consistent with their uncertainties. 

A simultanoueous use of A, B and C (line flagged with NSK+ABC) provokes an upward 
shift by 4 1CT 10 , as already observed in Tabled The behavior of the combination under the fit 
illustrates, once again, that the shape distortions exhibited by each of A, B and C are different 
and that the compensation performed by the fit degrades the description of the e + e data. This 
produces the increased value for a M , as also observed in the previous Subsection. 

21 And keeping the (X/rj) 2 term included in the % 2 expression, as stated above. 
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Data Set 


Fit Solution 


Statistical Information 




-inlO / \ 
10 U a ^{7171) 


X 2 /dof 


Probability 


iNofv (e e ) 


ooy.oi ± i.oz 


525.10/597 


98.4% 


NSK +A sh B sh C sh 


359.62 ± 1.51 


605.90/679 


97.9% 


NSK + A (Xt — 28-39) 


362.24 ± 1.52 


568.60/632 


96.6% 


NSK + B (x; = 32.59) 


360.68 ± 1.47 


558.88/614 


94.6% 


MC17 ■ /"< /,,2 if) |i\ 

JNbls. + L, (% T = jy.li) 


ooU.Oz ± 1.00 


565.24/624 


95.5% 


JNbK. + A r> 


oo4.4o ± l.o4 


648.60/680 


80.1 % 


NSK + KLOE 


359.22 ± 1.32 


652.92/657 


51.9% 


NSK + KLOE + A s/l B s/l C s/l 


358.52 ± 1.32 


741.50/739 


46.7% 


NSK + KLOE + A B C 


364.04 ± 1.25 


792.28/740 


8.9% 


NSK + KLOE + A(x 2 T = 49.59) 


361.55 ± 1.31 


708.90/692 


32.0% 


NSK + KLOE + B(x 2 T = 34.06) 


360.19 ± 1.16 


688.02/674 


34.6% 


NSK + KLOE + C ixl = 39.89) 


360.11 ± 1.32 


693.40/684 


39.3% 



Table 6: Contributions to 10 10 a /J (7T7r) from y/s G [0.630,0.958] GeV/c. r data set configu- 
rations are considered together with the e + e~ data and the appropriate subset of decay partial 
widths. The meaning of A, B, C, and A sft , B sh , C sh is given in the text. The contribution of 
single r subsets to the total x 2 is indicated when relevant by (x 2 - — ' ")• F° r lineshape fits 
without KLOE data, we have xl = 18-22 for A, xl = 28.35 for B and xl = 31.88 for C. 

4.2.3 Including KLOE Data 

Up to now, we have not introduced the KLOE lfi~8l data into the data sets submitted to a 
fit. Indeed, even if considered acceptable, its best x 2 / n points looks large [1J. It is nevertheless 
instructive to point out explicitly its effects. Table 2 in HI indicates that three among the 
rescaling coefficients can be cancelled out ; one only has to let vary those corresponding to the 
global scale (e ) and to the acceptance correction (e 2 ) lfl8l . 
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This provides the results given in the lower half of Tabled Even if the probability decreases 
due to the intrinsically high (minimum) x 2 of the KLOE data sample, one observes a good 
consistency of the a^nn) value derived from fitting the (rescaled) KLOE data together with 
all e + e~ Novosibirsk data. One also starts getting much improved uncertainties. One may note 
the good probability when using all NSK data together with KLOE and A sh B sh C sh and the 
improved uncertainties on a^irn). When using A, B and/or C , one finally observes the same 
trend as in the upper part of Tabled which lead us to be as cautious in this case. 

Nevertheless, one may conclude that KLOE data do not degrade the expectations from NSK 
data, whatever is done with r data. 

In order to assess the effects of a global fit, it is also interesting to compare the result for 
NSK+KLOE in Table [6] with the corresponding average published in 031 : 

GlobalFit : ^(tttt) = 359.22 ± 1.32 Average^ : a^inr) = 358.51 ± 2Al exp 

in units of 1CT 10 . As for NSK data alone, the central value from our fit is slightly larger than 
the estimate from [14J and the uncertainty is significatively improved. 

4.2.4 A Partial Summary 

In summary, the picture exhibited when using the r data samples looks somewhat intricated. 
We have shown above, especially in Subsection 13.11 that lineshape distortions and absolute 
scale of spectra are intimately related. 

However, the present analysis tends to indicate that an important part of the reported dis- 
crepancy between e + e _ and r based estimates of g — 2 is related with a more or less difficult 
way to account for isospin breaking effects differentiating the p ± and p° lineshapes. This con- 
clusion is enforced by comparing the behavior of the ALEPH data set with those of BELLE 
and CLEO : the distortions of the former set can be accounted for within our model, providing 
a reasonable accord with VMD expectations. 

As for BELLE and CLEO, genuine mass and width differences compared to p°, which suc- 
cessfully work with ALEPH data, do not avoid a residual distortion effect which has important 
consequences while estimating a M . 

Interestingly, the normalized r spectra do not reveal the same problem. This is certainly 
due to the fact that a free rescaling allows to disconnect lineshape distortions as parametrized 
by mass and width differences (Sg and 5m 2 in our model) from the absolute scale. This indi- 
cates that the distortions exhibited by BELLE and CLEO are not understood. Of course, one 
cannot exclude that the agreement with ALEPH is purely accidental and that a more refined IB 
mechanism might have to be considered. 

Figure [8] displays graphically the main conclusions discussed in this Section. The vertical 
line centered at the value derived by fitting all Novosibirsk data is drawn to guide the eye. 
Comparing our estimates with the experimental data reported there is also interesting by itself. 
This indicates that individual estimates of a M provided by each of the r data sets are close to 
the e + e~ based estimate. The combined fit of these, instead, exhibits a clear discrepancy, as 
shown by the data points flagged by NSK + A+B+C , NSK + Adm,d g +^dm,dg+Cdm,d g or with 
KLOE. Interestingly, the comparison of the various data points exhibits the same trend. 
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4.2.5 Comments On Our Model Uncertainties 

Uncertainties intrinsic to our model may have to be estimated. For this purpose, having 
alternative models able to cover the same scope than our extended model would be desirable. 
Unfortunately, such models do not seem to exist which could cope with detailed descriptions 
of several decay channels over our whole fitting range. Some interesting attempts have been 
initiated relying on Chiral Symmetry and the properties of the Roy equations 115611571 ; however, 
nothing final is presently available. 

In our case, one cannot make estimates by changing, for instance, the p parametrization 
(Gounaris-Sakurai, Breit-Wigner, Kuhn-Santamaria . . . ) as sometimes done elsewhere. In- 
deed, the form factor lineshape is intrinsic to the model, including the IB vector meson mixing. 
However, the difference between the central values of our estimates and the corresponding ex- 
perimental estimates lfT9ll20ll2"TTl may give some hint on the bound for our model uncertainty. 

4.3 Comparison with a M Estimates Collected by the ISR Method 

As a summary, our preferred final result is derived using all e + e~ data collected at Novosi- 
birsk : 

a^irir; 0.630 -=- 0.958 GeV) = [359.31 ± 1.62 exp ] 10~ 10 (22) 

and corresponds to an (underlying) fit probability of 98.4%. 

This result can be compared with the recent estimate published by KLOE ll54l : 

a^vrvr; 0.630 < m n7r < 0.958) = (356.7 ± 0.4 ± 3.1) 10~ 10 

based on a newly collected data set ifTBl [58l . This estimate may look a little bit low, but is 
consistent with ours at ~ 0.7 a level. 

On the other hand, the BaBar Collaboration has recently finalized a new data set [59], also 
collected using the ISR method. The contribution it provides to a^(mr) in the canonical s 
interval can be found in 11601 : 

a^nn; 0.630 + 0.958 GeV) = (365.2 ± 1.9 ± 1.9) 10" 10 = (365.2 ± 2.7 tot ) 10~ 10 

which supersedes a preliminary result presented at the Novosibirsk Conference TAU08[61| 
(a^iim) = (369 ± 0.8 ± 2.9) 10~ 10 ). Even if larger than Eq. (|22l) . this estimate is consistent 
with ours at the 1.9 a level. 

Comparing our result with these recent estimates clearly illustrates the advantage of intro- 
ducing information beyond the e + e _ — > n + ix~ cross section and a framework which encom- 
passes annihilation processes and partial width decays of light mesons. Indeed, the IB schemes 
needed can be calibrated consistently on the data set. The fact that the HLS framework does 
not presently go much beyond the <\> mass region is certainly a handicap, but does not prevent 
interesting improvements. 

4.4 Contributions to a fi Up To 1 GeV 

The discussion above leads us to consider the best motivated results for a M as derived from 
fit to the Novosibirsk e + e~ data samples taking into account the photon hadronic VP (see Eq. 
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(36) in[Q]|). A challenging choice might be to include the r normalized spectra, which have 
been shown to yield a satisfactory description simultaneously with all e + e~ data. We have 
seen that BELLE and CLEO seem to meet problems with their shape distortions which are 
not clearly identified. However, as Ag m; s g is reasonably well understood, one may consider its 
results for a M . 

Using these data samples, one can estimate several contributions^ to a M up to 1 GeV. They 
are given in Table U\ We thus provide the results obtained while fitting without any r sample 
in the first data column. The last data column, instead, shows the influence of KLOE data 
MM- 



Process 


NSK (no t) 


NSK + A B C 


NSK + Asmjg 


NSK + KLOE + A B C 


7T + 7T~ 


492.02 ± 2.24 


491.66 ± 1.98 


496.20 ± 1.85 


492.24 ± 1.79 


7T°7 


4.53 ± 0.04 


4.53 ± 0.04 


4.55 ±0.04 


4.51 ±0.05 


vn 


0.17 ±0.01 


0.17 ±0.01 


0.17 ±0.01 


0.17 ±0.01 


ii 


0.01 ±0.00 


0.01 ±0.00 


0.01 ±0.00 


0.01 ±0.00 


7r + 7r~7r° 


36.99 ±0.55 


36.97 ±0.56 


36.83 ±0.56 


37.07 ±0.55 


Total 


533.72 ± 2.32 


533.34 ±2.07 


537.76 ±2.08 


534.01 ±1.90 


Fit Probability 


98.4% 


98.0% 


95.3% 


50.8% 



Table 7: Contributions to 10 a M from thresholds up to 1 GeV/c for various processes fitting 
the data sets indicated on the first line ; by NSK we mean the set of all e + e~ data except for 
KLOE. The errors provided merge the reported statistical and systematic uncertainties. 



One can note that the estimate based only on e + e~ annihilation data is consistently im- 
proved by including within the fitted data sets the r lineshapes from ALEPH, BELLE and 
CLEO. Instead, using Ks m ,Sg produces a shift by 4 units, while reducing the uncertainty in the 
same way than the r lineshapes. Whether this estimate should be prefered, is an open question, 
awaiting a better understanding of the r spectrum distortions. On the other hand, including 
KLOE data does not provokes differences with using only the Novosibirk data, but instead 
provides improved uncertainties. 

The always (negligible) contribution of the annihilation process 7/7 is a prediction which is 
entirely determined by the 777 decay channel and our model which implies a tight correlation 
between the 777 and 7/7 final states ([TJ. So close to its threshold, this contribution could have 
been expected small ; it should become obviously larger while including the region. 

22 Actually, one could have produced, as well, these contributions up to about 1.05 GeV/c in order to include 
the cj> region. As our fitted region includes the <f> peak, our results would be as reliable. However, the largest (KK) 
contribution involving the <j> region is presently left aside because of the issue raised by [5] and still unsolved. 
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As a final remark, one may conclude that the global fit method indeed performs as expected. 
Moreover, as far as we know, our method is the first one which can associate a probability to 
the various reported contributions to the muon anomalous magnetic moment. Of course, the 
quoted probabilities refer to the fits underlying the estimates and not the estimates themselves, 
as these are only derived from computations using the fit results (parameter values and full 
error covariance matrix). One may infer that our prefered numerical results for a M - any of the 
first two data columns in Table [7J- should increase the discrepancy between the prediction and 
the BNL measurement of the muon g — 2. Improved estimates can be expected from using also 
the new data samples based on the ISR method collected by KLOE and BaBar. 

5 Summary And Concluding Remarks 

In the first part of this study [1], we defined the Extended HLS Model and reminded the 
mechanisms implementing U(3)/SU(3)/SU(2) symmetry breaking, noticeably the vector me- 
son mixing provided by breaking Isospin symmetry. This was shown to provide a satisfactory 
simultaneous description of all e + e~ data sets considered. The annihilation channels success- 
fully covered (n + n~, 7r°7, 777 and 7r°7r + 7r~) represent a large amount of data. It is also the 
largest set of annihilation channels simultaneously analyzed so far. The present work confirms 
that the dipion spectrum in r decay is also in the reach of this model. 

The r data sets collected by the ALEPH, BELLE and CLEO Collaborations have been care- 
fully examined in isolation and combined within the context of a global fit performed together 
with a large set of e + e~ annihilation data. For this purpose, we found appropriate to introduce 
additional Isospin breaking (IB) effects providing distortions of the lineshape compared to 
p°. This was partly considered in our former study JH, but totally ignored in preliminary ver- 
sions of the present work like ll62l . This additional mechanism (IB shape distortions) consists 
of a coupling difference 5g of the and p° to a pion pair and a mass difference Sm 2 between 
these two mesons. 

Each of the ALEPH (A), BELLE (B) and CLEO (C) experiments reports on a global scale 
uncertainty for their dipion spectra. This scale uncertainty is much smaller for A (0.51%) than 
for B (1.53%) or C (1.75%). The absolute scale of each r spectrum can then be modified 
by a factor 1 + Xa/b/c, including a fit parameter constrained, for each of A, B and C, by 
its scale uncertainty. However, it has been proved that, in minimization procedures, 5g, Sm 2 
and these scales are tightly correlated. Therefore, one cannot interpret Xa/b/c 7^ as purely 
reflecting biases on the measured branching ratios B(r — >■ n ± n°i i ) = B nw . To be explicit, 
the connection between pure lineshape distortion parameters and absolute scale, which is the 
subject of Subsection 13 . 1 1 prevents to tag reliably an undoubtful numerical value for some 
rescaling. 

It has been shown - see Table Q] - that the ALEPH data are in accord with VMD expec- 
tations, provided some shape distortions (5g, 5m 2 ) are implemented. In this case, = is 
even found quite acceptable. This means that, relying on ALEPH data only, the presently ac- 
cepted value for B nn is not contradicted by our analysis. This result is tightly connected with 
the fact that the lineshape distortions of the ALEPH spectrum is well accounted for by our 
parametrization (5g, Sm 2 ). 

The picture is not that clear with BELLE and CLEO, indicating that some residual problem 
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survives, mostly visible in these data sets ; this could well be due to having a too simple- 
minded shape distortion mechanism ; however, this could also reflect specific systematics in 
these experiments or, as well, a real physical problem revealed by their larger statistics. A 
rescaling of their B nn downwards is found, however, in reasonable accord with their reported 
uncertainties. 

Moreover, the lineshapes of the normalized spectra provided by A, B and C happen to 
exhibit no problem at all when fitted together with the largest possible set of e + e~ annihilation 
data. This is certainly due to having dropped out the correlation between absolute scale and 
lineshape distortions. In this case, one even finds no need to introduce significant IB shape 
distortions. The r lineshapes and the e + e~ are optimally fitted as clear from the individual x 2 
given in Table |3] 

Collecting all pieces of information in this study, one cannot consider the existence of a 
fundamental issue, at the pion form factor level, between e + e" and r data. Indeed, Table 
|3j which summarizes our simultaneous description of e + e~ annihilation and r decay data, 
displays a fit quality of about 80% with quite acceptable parameter values. This may imply as 
well that B W7r is slightly overestimated or that some additional systematics affect the 1 /NdN/ ds 
spectrum in some experiments. 

The most important aim of the present work was to study the effects of a global fit on the 
estimation of hadronic contributions to a M , the muon g — 2, from the various thesholds up to 
1 GeV. For this purpose, our treatment of statistical and systematic errors accounted as closely 
as possible for the information provided for each data sample. 

One first has checked that the fit performs as one could expect on individual e + e~ — > 
data samples. Comparing our fit results with individual experimental information [fT9l l20l |2TI 
about the contribution of the reference region (0.630, 0.958) GeV to a fl (n + n^) is already in- 
teresting : Table 0] shows the information returned by the fits is in good agreement with the cor- 
responding experimental information. Performing the global fit with these data samples, where 
systematics have certainly been considered with great care, provides quite reasonable central 
values for the individual estimates of a^ix+Tx^) and shrinked uncertainties ; this shrinking be- 
comes noticeable when fitting simultaneously all e + e~ data from |fT9ll20"ll2Tl : The uncertainty 
is reduced from 3.02 to 1.72 in units of 10~ 10 , i.e. better than a 40 % improvement. 

The next step was to include the older e + e~ — > data samples 11491 |5Tfl , where 

systematics are certainly not as well controlled as in [fP9l l20l l2Tfl . One thus gets a shift 
Aa At (7r + 7r~) ~ —1.9 (in units of 10~ 10 ) while leaving the uncertainty nearly unchanged. In- 
cluding also all the e + e _ — > (n° /rj)^ and e + e~ — > data samples does not produce 
significant modifications, showing that the effects of systematics when combining data sam- 
ples of various qualities may prevent improvements. However, the corresponding combined fit 
does not degrade the information ; it rather allows to check the stability of the estimate. The 
final uncertainty improvement is at the level of 46 %. 

At this point where all information on e + e" annihilations - except for KLOE ISR data - 
are considered, one ends up with : 

a^n + n-; 0.630 0.958) = 359.31 ± 1.62 

(in units of 10~ 10 ) where the error combines systematic and statistical uncertainties and ac- 
counts for the sample-to-sample correlations. The probability of the underlying fit to the data 
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is 98.4%. One has shown that these high probability values are a normal consequence of too 
conservative estimates of systematics into given data sets. This has been shown not to hide bad 
descriptions of some data subsets. 

Including KLOE data does not modify the picture, except for the probabilities which may 
become very low, reflecting the intrinsic large minimum \ 2 for this data set. 

Adding the various r data to the set of e + e _ data samples, is also worth being made step- 
wise, in order to substantiate problems and specific properties of each of the A (ALEPH), B 
(BELLE) and C (CLEO) data sets. We have first compared (see Table 0) our reconstructed val- 
ues for a M with the r based estimates of [fT4| for the reference region 0.630 4 0.958 GeV [|55l . 
The central values are found in agreement while the uncertainties are importantly shrinked. 
This is clearly an effect of our global fit/modelling. One should also note that the parameters 
provided by the simultaneous fit to A, B, C and the e + e~ data provides a value for a M larger 
than for each of A, B and C separately. Interestingly, this property is also exhibited by the 
combined experimental estimate proposed by the authors of IT141I551 . 

In this case, we show that introducing the A, B and C lineshapes inside our fit procedure 
allows to confirm the value for a /i (7r + 7r~; 0.630 4 0.958) derived using only e + e~ data and 
some more shrinking of its uncertainty : 



Both results correspond to good probabilities of the underlying form factor fits. 

However, our study leads us to conclude that some distortions of the lineshape, different 
for A and B/C, are at work which still need to be understood. Whether one is faced here 
with an incomplete IB distortion modelling, with unaccounted for systematics, or with some 
external physical effect, is an open issue. Until this issue is clarified, going much beyond the r 
lineshapes for a^(n + ix~) estimates looks hazardous. 

Our results concerning a M (7r + 7r~; 0.630 4 0.958) are summarized in Figure[8j This clearly 
illustrates that shape distortions and systematics are not completely understood, even if at the 
form factor level the picture is more optimistic ; this shows that a M estimates are a more sensi- 
tive probe to differences than fit probabilities of the underlying form factors. 

Our tendency for now is to prefer providing our final results concerning the various contri- 
butions to a M from thresholds to 1 GeV, considering all NSK e + e~ data samples together with 
r lineshape data. The KLOE data set does not seem to modify the picture unreasonably. This 
information is the matter of Table U\ 

One may also conclude that the global fit method performs as could be expected and is 
a useful tool in order to examine the consistency of various data sets covering various pro- 
cesses related by physics. As a tool, it also permits improved estimates of the various hadronic 
contributions to the muon g — 2. 

For this purpose, one should stress that better data, with better controlled systematics in all 
decay channels, and not only for the n + ix~ final state, may be valuable. Indeed, the physics 
correlations between the various final states in e + e _ annihilations may conspire (better than 
presently) with each other in order to provide improved values for all contributions. 



a M (7r + 7r - ; 0.630 4-0.958) 
Using additionally the KLOE data leads to : 

a M (7r + 7r~; 0.630 4 0.958) 



[359.62 ± 1.51] 10 



[358.52 ± 1.32] 10 



-10 
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Theoretical developments on the inclusion of scalar mesons and higher mass vector mesons 
within VMD-like frameworks are also desirable, as this could increase the underlying physics 
correlations between the various final states accessible in e + e~ annihilations. Finally, under- 
standing the issue raised by [5] about the coupling of the meson to KK pairs, may well be 
an important task in order to estimate reliably the region contribution to the muon g — 2. 
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Figure 1: Residual distributions of ALEPH Q3), BELLE El and CLEO (E) data sets in fits 
with mass and width differences between the charged and neutral p mesons. Upmost Figure 
includes fit rescaling factors contrained by the experimental uncertainty on B nn ; downmost 
Figure corresponds to assuming no rescaling. The fit region is bounded by 1 GeV/c. 
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Figure 2: Fit to the r spectra from ALEPH [15 J, BELLE H3 and CLEO data sets. The 
absolute normalization of the r spectra is completely free. 
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Figure 3: Residuals distribution in the fit to the r spectra from ALEPH lfT51 . BELLE [17J and 
CLEO [fT6ll data sets. The absolute normalization of the r spectra is completely free. 
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Figure 4: Residuals distribution in the fit to the r spectra from ALEPH lfT51 . BELLE [17J and 
CLEO lfT6l data sets. The residuals are normalized too the fit function values. The absolute 
normalization of the r spectra was left completely free. 
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Figure 5: Fit to the r spectra from ALEPH [15 J, BELLE H3 and CLEO OH data sets. The 
normalization is fit but constrained by the experimental uncertainties on B W7r (see text). 
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Figure 6: Residuals distribution in the fit to the r spectra from ALEPH |fl5ll . BELLE [17J and 
CLEO lfT6l data sets. The normalization is fit but constrained by the experimental uncertainties 
on B T7r (see text). 
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Figure 7: Residuals distribution in the fit to the r spectra from ALEPH lfT51 . BELLE [17J and 
CLEO lfT6l data sets. The residuals are normalized too the fit function values. The absolute 
normalization of the r spectra was constrained by the uncertainty on (see text). 
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Figure 8: Values of a^nn) x 10 10 integrated between 0.630 and 0.958 GeV. Data config- 
urations have been defined in Table @] for the first seven lines, and in Tables |5] and |6] for the 
following lines. The various notations used for the r data samples have been defined in the 
text. The points with dashed-dotted uncertainties are experimental values provided by the 
experiments lfT9l [201 12TH orinlfl4l. 
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